H3K4me1 recruits DNA repair proteins in plants.

DNA repair proteins can be recruited by their histone reader domains to specific epigenomic features, with consequences on intragenomic mutation rate variation. Here, we investigated H3K4me1-associated hypomutation in plants. We first examined two proteins which, in plants, contain Tudor histone reader domains: PRECOCIOUS DISSOCIATION OF SISTERS 5 (PDS5C), involved in homology-directed repair, and MUTS HOMOLOG 6 (MSH6), a mismatch repair protein. The MSH6 Tudor domain of Arabidopsis (Arabidopsis thaliana) binds to H3K4me1 as previously demonstrated for PDS5C, which localizes to H3K4me1-rich gene bodies and essential genes. Mutations revealed by ultradeep sequencing of wild-type and msh6 knockout lines in Arabidopsis show that functional MSH6 is critical for the reduced rate of single base substitution mutations in gene bodies and H3K4me1-rich regions. We explored the breadth of these mechanisms among plants by examining a large rice (Oryza sativa) mutation dataset. H3K4me1-associated hypomutation is conserved in rice as are the H3K4me1 binding residues of MSH6 and PDS5C Tudor domains. Recruitment of DNA repair proteins by H3K4me1 in plants reveals convergent, but distinct, epigenome-recruited DNA repair mechanisms from those well described in humans. The emergent model of H3K4me1-recruited repair in plants is consistent with evolutionary theory regarding mutation modifier systems and offers mechanistic insight into intragenomic mutation rate variation in plants.

Structural modeling of ion channels using AlphaFold2, RoseTTAFold2, and ESMFold.

Ion channels play key roles in human physiology and are important targets in drug discovery. The atomic-scale structures of ion channels provide invaluable insights into a fundamental understanding of the molecular mechanisms of channel gating and modulation. Recent breakthroughs in deep learning-based computational methods, such as AlphaFold, RoseTTAFold, and ESMFold have transformed research in protein structure prediction and design. We review the application of AlphaFold, RoseTTAFold, and ESMFold to structural modeling of ion channels using representative voltage-gated ion channels, including human voltage-gated sodium (NaV) channel - NaV1.8, human voltage-gated calcium (CaV) channel - CaV1.1, and human voltage-gated potassium (KV) channel - KV1.3. We compared AlphaFold, RoseTTAFold, and ESMFold structural models of NaV1.8, CaV1.1, and KV1.3 with corresponding cryo-EM structures to assess details of their similarities and differences. Our findings shed light on the strengths and limitations of the current state-of-the-art deep learning-based computational methods for modeling ion channel structures, offering valuable insights to guide their future applications for ion channel research.

Harnessing AlphaFold to reveal state secrets: Prediction of hERG closed and inactivated states.

To design safe, selective, and effective new therapies, there must be a deep understanding of the structure and function of the drug target. One of the most difficult problems to solve has been resolution of discrete conformational states of transmembrane ion channel proteins. An example is KV11.1 (hERG), comprising the primary cardiac repolarizing current, IKr. hERG is a notorious drug anti-target against which all promising drugs are screened to determine potential for arrhythmia. Drug interactions with the hERG inactivated state are linked to elevated arrhythmia risk, and drugs may become trapped during channel closure. However, the structural details of multiple conformational states have remained elusive. Here, we guided AlphaFold2 to predict plausible hERG inactivated and closed conformations, obtaining results consistent with myriad available experimental data. Drug docking simulations demonstrated hERG state-specific drug interactions aligning well with experimental results, revealing that most drugs bind more effectively in the inactivated state and are trapped in the closed state. Molecular dynamics simulations demonstrated ion conduction that aligned with earlier studies. Finally, we identified key molecular determinants of state transitions by analyzing interaction networks across closed, open, and inactivated states in agreement with earlier mutagenesis studies. Here, we demonstrate a readily generalizable application of AlphaFold2 as a novel method to predict discrete protein conformations and novel linkages from structure to function.

Toward Digital Twin Technology for Precision Pharmacology.

The authors demonstrate the feasibility of technological innovation for personalized medicine in the context of drug-induced arrhythmia. The authors use atomistic-scale structural models to predict rates of drug interaction with ion channels and make predictions of their effects in digital twins of induced pluripotent stem cell-derived cardiac myocytes. The authors construct a simplified multilayer, 1-dimensional ring model with sufficient path length to enable the prediction of arrhythmogenic dispersion of repolarization. Finally, the authors validate the computational pipeline prediction of drug effects with data and quantify drug-induced propensity to repolarization abnormalities in cardiac tissue. The technology is high throughput, computationally efficient, and low cost toward personalized pharmacologic prediction.

Elucidating molecular mechanisms of protoxin-II state-specific binding to the human NaV1.7 channel.

Human voltage-gated sodium (hNaV) channels are responsible for initiating and propagating action potentials in excitable cells, and mutations have been associated with numerous cardiac and neurological disorders. hNaV1.7 channels are expressed in peripheral neurons and are promising targets for pain therapy. The tarantula venom peptide protoxin-II (PTx2) has high selectivity for hNaV1.7 and is a valuable scaffold for designing novel therapeutics to treat pain. Here, we used computational modeling to study the molecular mechanisms of the state-dependent binding of PTx2 to hNaV1.7 voltage-sensing domains (VSDs). Using Rosetta structural modeling methods, we constructed atomistic models of the hNaV1.7 VSD II and IV in the activated and deactivated states with docked PTx2. We then performed microsecond-long all-atom molecular dynamics (MD) simulations of the systems in hydrated lipid bilayers. Our simulations revealed that PTx2 binds most favorably to the deactivated VSD II and activated VSD IV. These state-specific interactions are mediated primarily by PTx2's residues R22, K26, K27, K28, and W30 with VSD and the surrounding membrane lipids. Our work revealed important protein-protein and protein-lipid contacts that contribute to high-affinity state-dependent toxin interaction with the channel. The workflow presented will prove useful for designing novel peptides with improved selectivity and potency for more effective and safe treatment of pain.

Structural modeling of hERG channel-drug interactions using Rosetta.

The human ether-a-go-go-related gene (hERG) not only encodes a potassium-selective voltage-gated ion channel essential for normal electrical activity in the heart but is also a major drug anti-target. Genetic hERG mutations and blockage of the channel pore by drugs can cause long QT syndrome, which predisposes individuals to potentially deadly arrhythmias. However, not all hERG-blocking drugs are proarrhythmic, and their differential affinities to discrete channel conformational states have been suggested to contribute to arrhythmogenicity. We used Rosetta electron density refinement and homology modeling to build structural models of open-state hERG channel wild-type and mutant variants (Y652A, F656A, and Y652A/F656 A) and a closed-state wild-type channel based on cryo-electron microscopy structures of hERG and EAG1 channels. These models were used as protein targets for molecular docking of charged and neutral forms of amiodarone, nifekalant, dofetilide, d/l-sotalol, flecainide, and moxifloxacin. We selected these drugs based on their different arrhythmogenic potentials and abilities to facilitate hERG current. Our docking studies and clustering provided atomistic structural insights into state-dependent drug-channel interactions that play a key role in differentiating safe and harmful hERG blockers and can explain hERG channel facilitation through drug interactions with its open-state hydrophobic pockets.

A Novel Peptide Prevents Enterotoxin- and Inflammation-Induced Intestinal Fluid Secretion by Stimulating Sodium-Hydrogen Exchanger 3 Activity.

BACKGROUND & AIMS: Acute diarrheal diseases are the second most common cause of infant mortality in developing countries. This is contributed to by lack of effective drug therapy that shortens the duration or lessens the volume of diarrhea. The epithelial brush border sodium (Na+)/hydrogen (H+) exchanger 3 (NHE3) accounts for a major component of intestinal Na+ absorption and is inhibited in most diarrheas. Because increased intestinal Na+ absorption can rehydrate patients with diarrhea, NHE3 has been suggested as a potential druggable target for drug therapy for diarrhea. METHODS: A peptide (sodium-hydrogen exchanger 3 stimulatory peptide [N3SP]) was synthesized to mimic the part of the NHE3 C-terminus that forms a multiprotein complex that inhibits NHE3 activity. The effect of N3SP on NHE3 activity was evaluated in NHE3-transfected fibroblasts null for other plasma membrane NHEs, a human colon cancer cell line that models intestinal absorptive enterocytes (Caco-2/BBe), human enteroids, and mouse intestine in vitro and in vivo. N3SP was delivered into cells via a hydrophobic fluorescent maleimide or nanoparticles. RESULTS: N3SP uptake stimulated NHE3 activity at nmol/L concentrations under basal conditions and partially reversed the reduced NHE3 activity caused by elevated adenosine 3',5'-cyclic monophosphate, guanosine 3',5'-cyclic monophosphate, and Ca2+ in cell lines and in in vitro mouse intestine. N3SP also stimulated intestinal fluid absorption in the mouse small intestine in vivo and prevented cholera toxin-, Escherichia coli heat-stable enterotoxin-, and cluster of differentiation 3 inflammation-induced fluid secretion in a live mouse intestinal loop model. CONCLUSIONS: These findings suggest pharmacologic stimulation of NHE3 activity as an efficacious approach for the treatment of moderate/severe diarrheal diseases.

Energetic landscape of polycystin channel gating.

Members of the polycystin family (PKD2 and PKD2L1) of transient receptor potential (TRP) channels conduct Ca2+ and depolarizing monovalent cations. Variants in PKD2 cause autosomal dominant polycystic kidney disease (ADPKD) in humans, whereas loss of PKD2L1 expression causes seizure susceptibility in mice. Understanding structural and functional regulation of these channels will provide the basis for interpreting their molecular dysregulation in disease states. However, the complete structures of polycystins are unresolved, as are the conformational changes regulating their conductive states. To provide a holistic understanding of the polycystin gating cycle, we use computational prediction tools to model missing PKD2L1 structural motifs and evaluate more than 150 mutations in an unbiased mutagenic functional screen of the entire pore module. Our results provide an energetic landscape of the polycystin pore, which enumerates gating sensitive sites and interactions required for opening, inactivation, and subsequent desensitization. These findings identify the external pore helices and specific cross-domain interactions as critical structural regulators controlling the polycystin ion channel conductive and nonconductive states.

Opening of capsaicin receptor TRPV1 is stabilized equally by its four subunits.

Capsaicin receptor TRPV1 is a nociceptor for vanilloid molecules, such as capsaicin and resiniferatoxin (RTX). Even though cryo-EM structures of TRPV1 in complex with these molecules are available, how their binding energetically favors the open conformation is not known. Here, we report an approach to control the number of bound RTX molecules (0-4) in functional rat TRPV1. The approach allowed direct measurements of each of the intermediate open states under equilibrium conditions at both macroscopic and single-molecule levels. We found that RTX binding to each of the four subunits contributes virtually the same activation energy, which we estimated to be 1.70 to 1.86 kcal/mol and found to arise predominately from destabilizing the closed conformation. We further showed that sequential bindings of RTX increase open probability without altering single-channel conductance, confirming that there is likely a single open-pore conformation for TRPV1 activated by RTX.

Molecular determinants of µ-conotoxin KIIIA interaction with the human voltage-gated sodium channel NaV1.7.

The voltage-gated sodium (NaV) channel subtype NaV1.7 plays a critical role in pain signaling, making it an important drug target. Here we studied the molecular interactions between µ-Conotoxin KIIIA (KIIIA) and the human NaV1.7 channel (hNaV1.7). We developed a structural model of hNaV1.7 using Rosetta computational modeling and performed in silico docking of KIIIA using RosettaDock to predict residues forming specific pairwise contacts between KIIIA and hNaV1.7. We experimentally validated these contacts using mutant cycle analysis. Comparison between our KIIIA-hNaV1.7 model and the cryo-EM structure of KIIIA-hNaV1.2 revealed key similarities and differences between NaV channel subtypes with potential implications for the molecular mechanism of toxin block. The accuracy of our integrative approach, combining structural data with computational modeling, experimental validation, and molecular dynamics simulations, suggests that Rosetta structural predictions will be useful for rational design of novel biologics targeting specific NaV channels.

Elucidating Molecular Mechanisms of Protoxin-2 State-specific Binding to the Human NaV1.7 Channel.

Human voltage-gated sodium (hNaV) channels are responsible for initiating and propagating action potentials in excitable cells and mutations have been associated with numerous cardiac and neurological disorders. hNaV1.7 channels are expressed in peripheral neurons and are promising targets for pain therapy. The tarantula venom peptide protoxin-2 (PTx2) has high selectivity for hNaV1.7 and serves as a valuable scaffold to design novel therapeutics to treat pain. Here, we used computational modeling to study the molecular mechanisms of the state-dependent binding of PTx2 to hNaV1.7 voltage-sensing domains (VSDs). Using Rosetta structural modeling methods, we constructed atomistic models of the hNaV1.7 VSD II and IV in the activated and deactivated states with docked PTx2. We then performed microsecond-long all-atom molecular dynamics (MD) simulations of the systems in hydrated lipid bilayers. Our simulations revealed that PTx2 binds most favorably to the deactivated VSD II and activated VSD IV. These state-specific interactions are mediated primarily by PTx2's residues R22, K26, K27, K28, and W30 with VSD as well as the surrounding membrane lipids. Our work revealed important protein-protein and protein-lipid contacts that contribute to high-affinity state-dependent toxin interaction with the channel. The workflow presented will prove useful for designing novel peptides with improved selectivity and potency for more effective and safe treatment of pain.

Structural modeling of peptide toxin-ion channel interactions using RosettaDock.

Voltage-gated ion channels play essential physiological roles in action potential generation and propagation. Peptidic toxins from animal venoms target ion channels and provide useful scaffolds for the rational design of novel channel modulators with enhanced potency and subtype selectivity. Despite recent progress in obtaining experimental structures of peptide toxin-ion channel complexes, structural determination of peptide toxins bound to ion channels in physiologically important states remains challenging. Here we describe an application of RosettaDock approach to the structural modeling of peptide toxins interactions with ion channels. We tested this approach on 10 structures of peptide toxin-ion channel complexes and demonstrated that it can sample near-native structures in all tested cases. Our approach will be useful for improving the understanding of the molecular mechanism of natural peptide toxin modulation of ion channel gating and for the structural modeling of novel peptide-based ion channel modulators.

TRPV1 Opening is Stabilized Equally by Its Four Subunits.

Capsaicin receptor TRPV1 is a nociceptor for vanilloid molecules such as capsaicin and resiniferatoxin (RTX). Even though cryo-EM structures of TRPV1 in complex with these molecules are available, how their binding energetically favors the open conformation is not known. Here we report an approach to control the number of bound RTX molecules (0-to-4) in functional mouse TRPV1. The approach allowed direct measurements of each of the intermediate open states under equilibrium conditions at both macroscopic and single-molecule levels. We found that RTX binding to each of the four subunits contributes virtually the same activation energy, which we estimated to be 1.86 kcal/mol and found to arise predominately from destabilizing the closed conformation. We further showed that sequential bindings of RTX increase open probability without altering single-channel conductance, confirming that there is likely a single open-pore conformation for TRPV1 activated by RTX.

Elucidation of a dynamic interplay between a beta-2 adrenergic receptor, its agonist, and stimulatory G protein.

G protein-coupled receptors (GPCRs) represent the largest group of membrane receptors for transmembrane signal transduction. Ligand-induced activation of GPCRs triggers G protein activation followed by various signaling cascades. Understanding the structural and energetic determinants of ligand binding to GPCRs and GPCRs to G proteins is crucial to the design of pharmacological treatments targeting specific conformations of these proteins to precisely control their signaling properties. In this study, we focused on interactions of a prototypical GPCR, beta-2 adrenergic receptor (ß2AR), with its endogenous agonist, norepinephrine (NE), and the stimulatory G protein (Gs). Using molecular dynamics (MD) simulations, we demonstrated the stabilization of cationic NE, NE(+), binding to ß2AR by Gs protein recruitment, in line with experimental observations. We also captured the partial dissociation of the ligand from ß2AR and the conformational interconversions of Gs between closed and open conformations in the NE(+)-ß2AR-Gs ternary complex while it is still bound to the receptor. The variation of NE(+) binding poses was found to alter Gs α subunit (Gsα) conformational transitions. Our simulations showed that the interdomain movement and the stacking of Gsα α1 and α5 helices are significant for increasing the distance between the Gsα and ß2AR, which may indicate a partial dissociation of Gsα The distance increase commences when Gsα is predominantly in an open state and can be triggered by the intracellular loop 3 (ICL3) of ß2AR interacting with Gsα, causing conformational changes of the α5 helix. Our results help explain molecular mechanisms of ligand and GPCR-mediated modulation of G protein activation.

Adenylyl cyclase isoform 1 contributes to sinoatrial node automaticity via functional microdomains.

Sinoatrial node (SAN) cells are the heart's primary pacemaker. Their activity is tightly regulated by ß-adrenergic receptor (ß-AR) signaling. Adenylyl cyclase (AC) is a key enzyme in the ß-AR pathway that catalyzes the production of cAMP. There are current gaps in our knowledge regarding the dominant AC isoforms and the specific roles of Ca2+-activated ACs in the SAN. The current study tests the hypothesis that distinct AC isoforms are preferentially expressed in the SAN and compartmentalize within microdomains to orchestrate heart rate regulation during ß-AR signaling. In contrast to atrial and ventricular myocytes, SAN cells express a diverse repertoire of ACs, with ACI as the predominant Ca2+-activated isoform. Although ACI-KO (ACI-/-) mice exhibit normal cardiac systolic or diastolic function, they experience SAN dysfunction. Similarly, SAN-specific CRISPR/Cas9-mediated gene silencing of ACI results in sinus node dysfunction. Mechanistically, hyperpolarization-activated cyclic nucleotide-gated 4 (HCN4) channels form functional microdomains almost exclusively with ACI, while ryanodine receptor and L-type Ca2+ channels likely compartmentalize with ACI and other AC isoforms. In contrast, there were no significant differences in T-type Ca2+ and Na+ currents at baseline or after ß-AR stimulation between WT and ACI-/- SAN cells. Due to its central characteristic feature as a Ca2+-activated isoform, ACI plays a unique role in sustaining the rise of local cAMP and heart rates during ß-AR stimulation. The findings provide insights into the critical roles of the Ca2+-activated isoform of AC in sustaining SAN automaticity that is distinct from contractile cardiomyocytes.

Computational design of peptides to target NaV1.7 channel with high potency and selectivity for the treatment of pain.

The voltage-gated sodium NaV1.7 channel plays a key role as a mediator of action potential propagation in C-fiber nociceptors and is an established molecular target for pain therapy. ProTx-II is a potent and moderately selective peptide toxin from tarantula venom that inhibits human NaV1.7 activation. Here we used available structural and experimental data to guide Rosetta design of potent and selective ProTx-II-based peptide inhibitors of human NaV1.7 channels. Functional testing of designed peptides using electrophysiology identified the PTx2-3127 and PTx2-3258 peptides with IC50s of 7 nM and 4 nM for hNaV1.7 and more than 1000-fold selectivity over human NaV1.1, NaV1.3, NaV1.4, NaV1.5, NaV1.8, and NaV1.9 channels. PTx2-3127 inhibits NaV1.7 currents in mouse and human sensory neurons and shows efficacy in rat models of chronic and thermal pain when administered intrathecally. Rationally designed peptide inhibitors of human NaV1.7 channels have transformative potential to define a new class of biologics to treat pain.

Single amino acid residue mediates reciprocal specificity in two mosquito odorant receptors.

The southern house mosquito, Culex quinquefasciatus, utilizes two odorant receptors, CquiOR10 and CquiOR2, narrowly tuned to oviposition attractants and well conserved among mosquito species. They detect skatole and indole, respectively, with reciprocal specificity. We swapped the transmembrane (TM) domains of CquiOR10 and CquiOR2 and identified TM2 as a specificity determinant. With additional mutations, we showed that CquiOR10A73L behaved like CquiOR2. Conversely, CquiOR2L74A recapitulated CquiOR10 specificity. Next, we generated structural models of CquiOR10 and CquiOR10A73L using RoseTTAFold and AlphaFold and docked skatole and indole using RosettaLigand. These modeling studies suggested space-filling constraints around A73. Consistent with this hypothesis, CquiOR10 mutants with a bulkier residue (Ile, Val) were insensitive to skatole and indole, whereas CquiOR10A73G retained the specificity to skatole and showed a more robust response than the wildtype receptor CquiOR10. On the other hand, Leu to Gly mutation of the indole receptor CquiOR2 reverted the specificity to skatole. Lastly, CquiOR10A73L, CquiOR2, and CquiOR2L74I were insensitive to 3-ethylindole, whereas CquiOR2L74A and CquiOR2L74G gained activity. Additionally, CquiOR10A73G gave more robust responses to 3-ethylindole than CquiOR10. Thus, we suggest the specificity of these receptors is mediated by a single amino acid substitution, leading to finely tuned volumetric space to accommodate specific oviposition attractants.

Structural modeling of the hERG potassium channel and associated drug interactions.

The voltage-gated potassium channel, KV11.1, encoded by the human Ether-à-go-go-Related Gene (hERG), is expressed in cardiac myocytes, where it is crucial for the membrane repolarization of the action potential. Gating of the hERG channel is characterized by rapid, voltage-dependent, C-type inactivation, which blocks ion conduction and is suggested to involve constriction of the selectivity filter. Mutations S620T and S641A/T within the selectivity filter region of hERG have been shown to alter the voltage dependence of channel inactivation. Because hERG channel blockade is implicated in drug-induced arrhythmias associated with both the open and inactivated states, we used Rosetta to simulate the effects of hERG S620T and S641A/T mutations to elucidate conformational changes associated with hERG channel inactivation and differences in drug binding between the two states. Rosetta modeling of the S641A fast-inactivating mutation revealed a lateral shift of the F627 side chain in the selectivity filter into the central channel axis along the ion conduction pathway and the formation of four lateral fenestrations in the pore. Rosetta modeling of the non-inactivating mutations S620T and S641T suggested a potential molecular mechanism preventing F627 side chain from shifting into the ion conduction pathway during the proposed inactivation process. Furthermore, we used Rosetta docking to explore the binding mechanism of highly selective and potent hERG blockers - dofetilide, terfenadine, and E4031. Our structural modeling correlates well with much, but not all, existing experimental evidence involving interactions of hERG blockers with key residues in hERG pore and reveals potential molecular mechanisms of ligand interactions with hERG in an inactivated state.

Common Structural Pattern for Flecainide Binding in Atrial-Selective Kv1.5 and Nav1.5 Channels: A Computational Approach.

Atrial fibrillation (AF) is the most common cardiac arrhythmia. Its treatment includes antiarrhythmic drugs (AADs) to modulate the function of cardiac ion channels. However, AADs have been limited by proarrhythmic effects, non-cardiovascular toxicities as well as often modest antiarrhythmic efficacy. Theoretical models showed that a combined blockade of Nav1.5 (and its current, INa) and Kv1.5 (and its current, IKur) ion channels yield a synergistic anti-arrhythmic effect without alterations in ventricles. We focused on Kv1.5 and Nav1.5 to search for structural similarities in their binding site (BS) for flecainide (a common blocker and widely prescribed AAD) as a first step for prospective rational multi-target directed ligand (MTDL) design strategies. We present a computational workflow for a flecainide BS comparison in a flecainide-Kv1.5 docking model and a solved structure of the flecainide-Nav1.5 complex. The workflow includes docking, molecular dynamics, BS characterization and pattern matching. We identified a common structural pattern in flecainide BS for these channels. The latter belongs to the central cavity and consists of a hydrophobic patch and a polar region, involving residues from the S6 helix and P-loop. Since the rational MTDL design for AF is still incipient, our findings could advance multi-target atrial-selective strategies for AF treatment.

Structural and Functional Characterization of a Novel Scorpion Toxin that Inhibits NaV1.8 via Interactions With the DI Voltage Sensor and DII Pore Module.

Voltage-gated sodium channel NaV1.8 regulates transmission of pain signals to the brain. While NaV1.8 has the potential to serve as a drug target, the molecular mechanisms that shape NaV1.8 gating are not completely understood, particularly mechanisms that couple activation to inactivation. Interactions between toxin producing animals and their predators provide a novel approach for investigating NaV structure-function relationships. Arizona bark scorpions produce Na+ channel toxins that initiate pain signaling. However, in predatory grasshopper mice, toxins inhibit NaV1.8 currents and block pain signals. A screen of synthetic peptide toxins predicted from bark scorpion venom showed that peptide NaTx36 inhibited Na+ current recorded from a recombinant grasshopper mouse NaV1.8 channel (OtNaV1.8). Toxin NaTx36 hyperpolarized OtNaV1.8 activation, steady-state fast inactivation, and slow inactivation. Mutagenesis revealed that the first gating charge in the domain I (DI) S4 voltage sensor and an acidic amino acid (E) in the DII SS2 - S6 pore loop are critical for the inhibitory effects of NaTx36. Computational modeling showed that a DI S1 - S2 asparagine (N) stabilizes the NaTx36 - OtNaV1.8 complex while residues in the DI S3 - S4 linker and S4 voltage sensor form electrostatic interactions that allow a toxin glutamine (Q) to contact the first S4 gating charge. Surprisingly, the models predicted that NaTx36 contacts amino acids in the DII S5 - SS1 pore loop instead of the SS2 - S6 loop; the DII SS2 - S6 loop motif (QVSE) alters the conformation of the DII S5 - SS1 pore loop, enhancing allosteric interactions between toxin and the DII S5 - SS1 pore loop. Few toxins have been identified that modify NaV1.8 gating. Moreover, few toxins have been described that modify sodium channel gating via the DI S4 voltage sensor. Thus, NaTx36 and OtNaV1.8 provide tools for investigating the structure-activity relationship between channel activation and inactivation gating, and the connection to alternative pain phenotypes.